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Based on the result of an unzipping phase transition by a force in a model of DNA, we hypothesize that the 
DnaA-type proteins act as a pulling agent with a force slightly less than the critical force for unzipping. The dynamic 
(space-time) correlation of unzipping then drives the subsequent events of replication. In such a correlation driven 
scenario, there is no need of a replisome as a structural unit and this elusive replisome may not exist at all. 



I. INTRODUCTION 

A double stranded DNA (dsDNA) needs to be opened up for its replication ffl. Once done, the subsequent processes 
take place sequentially, like the binding of, say, helicases, single-strand-binding proteins etc, then the elongation of 
the new chain by polymerases etc. 

Recently it has been proposed |2j that the action of dnaA-type proteins could be thought of as a force acting at 
the initiation site (called "origin" ) . ( The function of DnaA is not to be confused with the "locomotive" or "sewing 
machine" action of the hypothetical replisome.) A simple model showed that there is a phase transition at a critical 
value of the force so that for forces less than the critical strength, DNA is a double-stranded object while for forces 
exceeding it, the chains get opened up. 

In this paper we give a short description of this model, the principal results and the subsequent developments. We 
point out the features that seem to have led to "confusions" . We also show how this simple model can be extended 
to incorporate some other features. 

The connection to in vivo replication and laboratory experiments are discussed in section III, where we propose 
that the correlation of unzipping is at the heart of the highly correlated collective process that goes by the name of 
replication. 

II. STATISTICAL MECHANICAL MODEL 
A. Model and results 

The model proposed originally is that of a pair of flexible polymers bound together by a short-range interaction 
and pulled at one end by a force We start with a homo-DNA, i.e. all base pairs identical. One end of the two 
strands are anchored or tied together and the other end is pulled by a force g. The Hamiltonian is given by 
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where rj(s) denotes the ci-dimensional position of a monomer on chain i (of elastic constant ei) at a contour- length s 
measured from the tied end (s = 0), r = ri — r 2 is the relative separation of the two chains at the same monomer s 
and V(r) is a short-range potential simulating the interaction of the base pairs of two strands. Any realistic potential 
can be chosen here. But our interest is in the effect of the force on the bound double-strand case, and so for analytical 
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simplicity we choose a delta function (or contact) potential, to be interpreted as a limit of a narrow square well. The 
strength of the well is chosen to be such that there is only one bound state in the problem. 

As usual in statistical mechanics, there are two possible ensembles, namely a constant force (g constant) ensemble 
and a constant end-point separation(r(./V)=constant) ensemble. For studies of phase transitions, it is useful to work 
in an ensemble of fixed intensive variable. We therefore choose the fixed force ensemble. 

For the zero force case, there is a critical unbinding transition at v — v c . The pulling force would like to align the 
strands in the direction of the force while the interaction would like to keep the strands together. There is a critical 
force below which the DNA remains in the dsDNA phase while for forces exceeding this critical value it gets unzipped. 

Quantitatively, the partition function can be evaluated by a transfer matrix approach with the contour length as 
the transfer direction. The free energy per unit length is given by the largest eigenvalue of the transfer matrix and a 
phase transition takes place whenever there is a degeneracy of the largest eigenvalue. 

The transfer matrix can conveniently be written in the form of a quantum Hamiltonian, albeit non-hermitian 
for a particle of co-ordinate r (CM behaviour is like a free chain or a free particle) 

H q ( g ) = ±(p + t g) 2 + V(r), (2) 

in units of h{= fceT) — 1 and mass = 1, with p as momentum. For long chains N — > oo the free energy is the 
ground-state energy of this non-hermitian Hamiltonian. A phase transition takes place whenever the ground state 
is degenerate. The analysis done in Ref. || shows that if the ground state energy ( i.e. the binding energy of the 
dsDNA per unit length or per base pair ) is E , then the critical force is given by 

9c = 2tM~ \v-v c \ l '\ 2 ' d \. (3) 

where the w-dependences [Q of Eq close to v c , for general d, is used. In absence of any force, DNA can be denatured 
either by changing temperature or by changing pH of the solution. In this respect the effective interaction parameter 
v is a better variable than the temperature itself. The generic phase diagram is shown in Fig. 1. 




FIG. 1. The g — v phase diagram. Thermal or pH driven denaturation without any force is the g — line. The DnaA 
protein is hypothesized to work in the shaded region just below the critical line. 

The degeneracy of the eigenvalues leads to a diverging length scale that could be associated, for g — > g c — . as the 
characteristic length of the opened Y-fork (the replication fork). This length scale measured along the strand of the 
DNA gives a measure of the number of monomers opened by the force. Denoting this number by m, we get 

(m) ~| g c - g \~" m with v m = 1. (4) 

(y m was denoted by v T in Ref Q].) 

B. Questions on the model and the results 

It is known for centuries that a certain minimum pulling force is needed to move an object on a surface with 
friction, or for adhesion, and so, is the occurrence of a critical force an obvious result? An understanding of this 
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threshold phenomena beyond phenomenological level is still in its infancy. Our results are to be seen in this context, 
and the model is different from the stick-slip models used in friction or adhesion problems. In fact, previous force 
measurements |] for unzipping of DNA have been interpreted in terms of equilibrium statistical mechanics, and not 
as a friction problem. 

A classical diatomic molecule or a bond (or a particle in a potential well) under a force does not show any phase 
transition. At a quantum mechanical level, a bond with a stretching force is still described by a hermitian hamiltonian 
and therefore will not show the phase transition of the previous subsection. It is the polymeric property or connectivity 
that is crucial for the unzipping transition. 

The result of Eq. ^ is not to be confused with the pulling of a polymer discussed in text books. In the case of a 
polymer pulled at the ends by a stretching force, there would be an extension of the chain and for a small force, the 
overall dimension of the polymer would scale linearly with the stretching force In contrast, we see that, in the 
dsDNA situation, a small force will have no major effect and the law in Eq. ^ is a dependence on the deviation from 
the critical force. 

Is continuum model a goo d approach? The answer to this question lies in the recognition of various length scales. 
The Hamiltonian of Eq. [la| can be discretized and studied with appropriate bond lengths and realistic V(r). We 
adopted a coarse-grained approach to focus on the phase transition behaviour. In other words, a continuum approach 
by itself is not a necessity to study the phase transition and other details of unzipping. 

Previous studies seem to indicate that DNA is better represented by worm-like chains but again a renormalization 
group argument shows that in the long length limit such a Hamiltonian will generate the elastic term |?j and therefore 
Eq. |la| is a representative of the universality class. In the same spirit the delta-function potential can be replaced 
by any realistic potential. But for the binding-unbinding ( melting) transition, it is only the integral of the potential 
that matters and so a square well or delta function is a valid starting point. 

Is the restriction to homo-DNA any good? It is worth remembering that the genetic code was deciphered from 
homoDNA's. For a controlled experiment to detect the unzipping transition, large homo-DNA's should be studied. 
The model is easily generalizable (see below) to consider specific base-sequence 



C. Subsequent developments 

After the original proposal of model for unzipping in Ref. several studies have been made. The dynamics of 
pulling a polymer from a potential well, i.e., the dynamics of unzipping in one dimension has been studied at a mean 
field level by Sebastian ||]. In addition to recovering the results of previous section, Ref |8) shows that the unzipping 
transition is also characterized by a diverging time scale. Denoting the characteristic time scale by t, one finds 

t~ISc-#r i/ ' witn v t = 1- (5) 

A divergence of time scales indicates a long-range correlation in time in dynamics of unzipping. 

DNA generally consists of inhomogeneous sequence of base pairs. Such a case in this model can be considered by 
taking the interaction energy to be monomer position dependent. The interaction term in Eq. [la| can be written as 
J ds b(s)v S(r), with 6(s) a variable depending only on monomer (base) position s describing the specific details 
of the sequence. A case of random b(s), a random interaction model (RANI model) has been considered earlier || 
for the zero force case where the randomness has been found to be marginal. In the case of unzipping transition, the 
randomness does not destroy the sharpness of the transition but changes the exponent v m defined in Eq. || to 
v m = 2. The dynamics is yet to be studied. 

The nature of the transition has also been considered. Maritan, Orlandini and Seno [|llj showed by Monte Carlo 
simulations and analytical methods that the noncrossing constraint on the chains can lead to a first order transition, 
and more interestingly to a "cold denaturation" . It has recently been argued [|l2| from a direct evaluation of the 
partition function that the transition for the Hamiltonian of Eq. [la] for d = 1 could be first order. It is known that 
depending on the reunion exponent, the phase transition for polymers could be first order Jl3[-p6|, and in the quantum 
context, the order of transition (Eq. ^) is determined by the normalizability of the "critical" wavefunction |]l7| , |l8| . 
Several exact results for the unzipping transition in various dimensions have been obtained in Ref. [Ell. 



D. Extensions 



The simple model does not take into account the self avoidance of the chains. Effects of self-avoidance has been 
studied in Ref. In a laboratory experiment, one can take Eq. [la| as a model for DNA in a theta-solvent §. In a 
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good solvent, the self and mutual avoidance of each chain can be introduced by an imaginary random potential and 
averaging the partition function J2(| . Adding a term to Eq. [la] 



imag 



, , ds [zV(r 1 ( S )) + zV(r 2 ( S ))] I (6) 

with V(r) as an annealed Gaussian random potential with zero mean and variance (V(r) V(r')} = uS(r — r'), the 
averaged partition function Z = J VIZ cxp(— H), with ovcrlinc denoting averaging over V, leads to a hamiltonian of 
the form of Eq. [la| with additional terms, H[ nt 



H int = dsds' V u*(r i («)-r i (fl')). (7) 



E 

i,i=l,2 

This represents excluded volume interaction for the chains. This imaginary potential of Eq. ^| would lead to a complex 
scalar potential in the quantum hamiltonian of Eq. |^. 

We have ignored the double helical nature of DNA. In actual replication, topoisomerases act when there is super- 
coiling. Once the chains are opened up in a region by DnaA forming a Y-fork, the remaining chain is brought back 
to the native sate by the topoisomerase. This gives a justification for the absence of the helical configuration in the 
hamiltonian in Eq. [la|. The topological linking number of the two chains can be described by the Gauss integral 
which is known to act like a real vector potential [ pO[ . It would be appropriate to take into account the linking as an 
independent quantity only if there is a change in the coiling in the process [ ^l| . 

A toy model is a quantum Hamiltonian corresponding to a polymer with winding as 



H q =-{p + A) 2 + V(v)+iV(r), (8) 



where A — Ao(r) + ig, with the real vector potential describing windings and the imaginary part, as before, the 
exerted force. Studies of the most general quantum Hamiltonian with both complex scalar and vector potential are 
needed. 



III. IMPLICATION FOR BIOLOGICAL REPLICATION 



We explain the hypothesis on the role of unzipping in replication. 

The replication process starts with a DnaA protein attaching itself, in interaction with the membrane, at the "origin" 
to start the Y-fork ( or "eye"). The next step is the binding of the various enzymes/proteins like single-strand binding 
(SSB) protein, helicase, topoisomerase, polymerases etc. Unlike the latter proteins/enzymes, the functionality of 
DnaA protein is not well understood p3 | . It is apparent that the whole process requires a strong correlation in space 
and time, whose source or origin is not obvious. In order to explain such correlations in subsequent processes, a 
structural unit "replisome" has been postulated which has never been isolated. 

We make the hypothesis that 

• The function of DnaA is actually to exert a pulling force (let us take this time t = 0) with a force which is close 
to, but slightly below, the critical force for unzipping. This force tends to form a Y-fork or an eye-type bubble 
(depending on the location of the origin). 

• The large length and time scales for the unzipping process, leads to a space time correlation in the initial 
nonequilibrium dynamics of unzipping. We postulate that this correlation controls the subsequent processes 
especially the dynamics of bindings of the subsequent enzymes/proteins. 

The correlation in the unzipping dynamics can be characterized quantitatively by the early time behaviour of, say, 
(m(t)m(i' — > 0)) where m(t) is the number of monomers unzipped at time t, or a more microscopic correlation 
(ri(s,t)rj(s' — > N,t' — > 0)), i,j being the spatial components. 
The dynamics of unzipping is described by a Langevin equation 

dr(s,t) SH 

where t represents time, H is given by Eq. [la], r\ is the thermal noise related to friction coefficient /i by 
(r}i(s,t)rjj(s' ,t')) = 2dijfj,kBT5(s — s') 5(t — t'). The binding of the next set of proteins in the opened region may be 
represented symbolically by a kinetic equation of the type 
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-= -C(s,t)f{c}, 



(10) 



where c denotes the concentration of the free enzymes in the bath, f{c} describes the dynamics in absence of any 
correlations, and C(s,t) stands for the initial or early time behaviour of the unzipping correlations. This is equivalent 
to making the rate constant dependent on the initial correlation p^] . This way of viewing the replication process 
as coupled system driven by the long-range space-time correlated unzipping dynamics is different from the prevalent 
contact-based approach. 

It is necessary to (a) study long homo-DNA's under well-controlled pulling force to map out the phase diagram 
^[H], and (b) probe jointly the correlation and binding kinetics in presence of a pulling force. 

A few appealing features of such a scenario are to be noted. Since the critical force is dependent on the structural 
details of the DNA, which is different for different species, our hypothesis naturally requires species variations of 
the DnaA-type protein for proper functionality. This protein should be such that it exerts a force slightly below 
the critical value for that particular DNA. Finally, in the proposed scenario, the dynamic correlation plays the most 
important role, and therefore a composite structure is not a necessity. In other words, the so-far-elusive hypothesized 
replisome need not exist as a structural unit. 
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